Robustness of multi-qubit entanglement in the independent decoherence model 
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We study the robustness of the GHZ (or "cat") class of multi-partite states under decoherence. 
The noise model is described by a general completely positive map for qubits independently coupled 
to the environment. In particular, the robustness of A^-party entanglement is studied in the large 
limit when (a) the number of spatially separated subsystems is fixed but the size of each subsystem 
becomes large (b) the size of the subsystems is fixed while their number becomes arbitrarily large. 
We obtain conditions for entanglement in these two cases. Among our other results, we show that 
the parity of an entangled state (i.e., whether it contains an even or odd number of qubits) can lead 
to qualitatively different robustness of entanglement under certain conditions. 
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I. INTRODUCTION 

Entanglement is believed to be a crucial resource for 
quantum information processing, e^., quantum comput- 
ing and quantum communication It is expected that 
practical realizations of quantum protocols will involve 
independent or collective manipulation of large-scale en- 
tanglement, i.e., entanglement distributed among N par- 
ticles, where N can be arbitrarily large. Unfortunately, in 
the absence of active intervention, entanglement is noto- 
riously susceptible to decoherence^^ ■ Intervention in the 
form of distillation protocols 0, Q^and error correcting 
codes is capable of restoring robustness to entangle- 
ment in the presence of decoherence. However, the effi- 
ciency of such methods strongly depends on the a priori 
robustness of the entangled state in question. This raises 
the natural question of the inherent robustness of entan- 
glement. Recent work has shown that for certain types of 
entangled states 0, 0, 111 or noisy preparation procedures 
there is indeed an (unexpected) inherent robustness. 
Here we continue this line of investigation and consider 
the inherent robustness of multi-qubit entangled states 
under a rather general model of uncorrelated decoher- 
ence. 

Let us here note that the notion of robustness of en- 
t an gled states has been used in a different context before 
ESIH in ■ However, in this paper by robustness we 
simply mean the ability of an entangled state to remain 
entangled in presence of decoherence. 

Besides the obvious practical importance of studying 
the effect of decoherence on multi-qubit entangled states, 
there is a fundamental interest as well. Entanglement be- 
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ing a microscopic property, one may ask how often macro- 
scopic entanglement is realized in the physical world? In 
other words, when the number of particles sharing an 
entangled state becomes very large, entanglement truly 
becomes a macroscopic property of the system itself. At 
the same time entanglement could become exponentially 
fragile, in the sense that an arbitrary small amount of 
noise can destroy the complete coherence between the 
superposed states. Yet, it is known that the set of sep- 
arable states is much smaller than the set of inseparable 
states ,23], which suggests that entanglement should be 
relatively common. These conflicting intuitions suggest 
that the question of the robustness of entanglement is a 
subtle one. 

The structure of A^-party entanglement is considerably 
more complex than a simple bipartite scenario where en- 
tanglement is only distributed among two subsystems. 
For recent results on multipartite entanglement and it's 
measures one can see 0|. An important no- 

tion is the partitioning of the system into 2 < M < N 
parties {M -partitioning) , where each of the M parties of 
several particles is considered to be a single system with 
a higher dimensional Hilbert space. A related notion 
is M-distillability: some M -party pure entanglement can 
be obtained by local operations and classical communica- 
tions (LOCC). A necessary condition for M-distillability 
is that the M -partitioned state is non-positive under par- 
tial transposition (NPPT) across all bipartite cuts |l7j . 
Thus to obtain information about distillability of an A'^- 
party state it suffices to study entanglement properties 
across all possible bipartite cuts. It is important to real- 
ize that if there is distillable entanglement between every 
pair then the whole state is M -distillable as any multi- 
partite entangled state can be prepared given sufficient 
bipartite entanglement between every pair. Note, that 
the original A^-party state not being distillable when all 
A^ parties are separated, does not rule out the state being 
distillable for some M-partitioning. 

To address the issue whether entanglement can also 
be viewed as a macroscopic property it is necessary to 
study the limit of large A'. If we assume a democratic 
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partitioning for simplicity, i.e., that N parties are divided 
into M groups such that Mk = N where k is the num- 
ber of particles in each group, then a natural question 
is, out of all possible Af-partitionings, which partition 
exhibits maximal robustness? Does entanglement exist 
when N ^ oo7 How does large scale entanglement be- 
have, for instance, when the size of each partition be- 
comes macroscopic while keeping the number of parti- 
tionings fixed? Such questions have been studied where 
the noise model was described by a depolarizing channel 
H, . We address these questions in the context of a 
rather general decoherence model in this work. 

The structure of the paper is as follows. In Sec. ITU we 
introduce the decoherence model: a general completely 
positive map under the assumption that each qubit is 
independently coupled to the environment. In Sec. lIIII we 
preview, for convenience, the main results of this work. 
Section IIVI finds the resulting density matrix under the 
action of decoherence. Section deals with some of the 
useful properties of the noisy density matrix. In Sec. IVII 
we obtain the entanglement conditions that quantify the 
robustness of the output state. Section IVIII concludes 
and lists some open problems. 



II. DECOHERENCE MODEL 

In Refs. 0, 0] properties of large-scale entanglement 
in the presence of the depolarizing channel were stud- 
ied. Ref. studied TV-particle GHZ (also known as 
"cat") states and compared them to W-states and 
spin-squeezed states. Ref. f?| studied the rather general 
class of graph states, which includes GHZ and cluster 
states. In Ref. the robustness of symmetric entangled 
states subject to particle loss was studied. Here we focus 
on A^-particle cat states, but considerably generalize the 
decoherence model. Results corresponding to the widely 
used depolarizing channel or dephasing channel can be 
reproduced as special cases of our model. 

An iV-qubit cat state is of the form I'l')^ = 

1/V2(|0)^^ + |1)®^), where |0),|1) are the +1,-1 

eigenstates of the Pauli az matrix. We assume that 
qubits are individually coupled to the environment. The 
action of the noisy channel on the mth qubit is described 
by a completely positive map of the form: 

3 

P — ' P ^ S,n (p) = TTop + ^ 7r,cr™pcr™ (1) 

1=1 

where the Tr{t) are probabilities (tt^ > 0, X]^=o ""i — l)i 
and the cr^ arc the Pauli matrices (cti = a^, etc.). The 
depolarizing channel and dephasing channel are the spe- 
cial cases TTi = (1 — 7ro)/3 (i = 1, 2, 3), and m = n2 ~ 0, 
respectively. By studying this rather general model we 
hope to provide more physical insight into the factors af- 
fecting the entanglement of a multi-qubit system in pres- 
ence of decoherence. 



III. SUMMARY OF MAIN RESULTS 

Before launching into the analysis, let us give a brief 
preview of the results obtained. We quantify the entan- 
glement robustness of iV-particle cat states by establish- 
ing sufficient conditions for the state to be entangled. 
Specifically, we obtain the conditions for M-distillability, 
for all 2 < M < N. These conditions are obtained for 
both finite and infinite N. 

• The robustness of the entangled state depends on 
sums and differences of the probabilities n and not 
on the actual probabilities. 

• The parity of the entangled state, i.e., whether the 
state is composed of an even or odd number of 
qubits, can lead to different qualitative behavior 
under certain conditions. 

• Macroscopic entanglement is found to be more ro- 
bust when distributed among higher dimensional 
subsystems, while keeping the number of spatially 
separated parties unchanged. 

• The most robust partition is found to be the bi- 
partite one where we have equal (TV is even) or 
approximately equal (TV is odd) number of qubits 
on both sides. 



IV. EXPRESSION FOR FINAL STATE 

We begin by studying the action of the noisy channel 
on the input state |^')jy ■ 

Pn ^ Pn = SiS2---Sn (pn) (2) 
The input pure state pN — |^')7v (^1 can be written as 

PN = l[m (oi)®^ + (ii>(ir^ 
+ (io)(ii)«^ + (ii) (oi)n 

^livT + vT+vT+vfo''] (3) 
where \i) {j\ = riv,hj 0,1. 

The action of the noise channel on the operators 
making up the density matrix is 

S{r]oo) = arjoo + &?7ii = Aqq, 
5(7711) = 6?7oo + a?7ii = All, 
S'(77oi) = crioi + driiQ = Aqi, 
S{r]io) ^ cr]io + drioi = Aio, (4) 

where 

a = TTo -I- TTs , 6 = TTl + 7r2 , C = TTo - TTs , d = TTl - 7r2 (5) 
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The coefficients a, b are positive wliereas tfie coefficients 
c, d can be eitlier positive or negative. Note that the 
grouping of ttq with tt^ and of tti with 7r2 is a result of 
our choice to work in the |0), |1) basis. The final density 
matrix can now be written as 

Pn = \ [AoT + Ar + A„T + Ar ] (6) 
Expanding the operator terms one obtains, 

N-l 



(El 



'AT 



+ permutations] 



(7) 



where 



denotes a ket containing exactly j 



Os and the suffix p denotes permutation. For a given j 
there are (^) permutations. 

Let us now define the A^-qubit "cat-basis" , consisting 
of 2^ states. First, 

|vl/±) = i=(|00...0)±|ll...l)). 

Other than these two states, we will call any other state 
a member of the fc-group if it is of the form: 



0;(7V- fc) : 1) ± |fc : 0; (TV - A:) : l)) 



where a = f (a,b,c,d). The index k stands for the 
group and the index j corresponds to the different states 
due to permutation of indices that belong to the same 
group. Comparing the two representations Eq. |(7J) and 
Eq. @ of the final density matrix, and noting the fact 



that ( *+ 



kj 



P 



a^,^ remain 



unchanged under local depolarizing, one can show that 
a^. +a'i^j = {a^b^~^ + b^a^'^) /2 which is independent 
of j. 



V. BASIS CONDITION FOR ENTANGLEMENT 
AND SOME USEFUL LEMMAS 

For any arbitrary A^-qubit density matrix p it has been 
shown that the state p is entangled if at least for one pair 
(fc, j), the following condition is satisfied 0, 



K*J|p|*J)-(*o~|p|*o")| > 



(10) 



kj 



Accordingly we compute the above quantities for our 
state. First observe that. 



an 



(*([|H*^)-(*o|/'|*(7) 

i[{(oo...oip;vioo...o) + (ii...i|p;viii-i>} 



+ {(00...0|pAr|ll.. 

-{(oo...o|p;vloo.. 

+{(oo...oip;viii- 

= 2Re(00...0|p^ 111 



l) + (ll...l|p^|00...0)} 

0) + (ll...l|p;v|ll-l>} 

1) + (ll...l|p;v|00...0)}] 
..1) (11) 



One can readily evaluate from the expression of the final 
density matrix the value of — a'^ , 



(8) 

This notation means that the first term of the super- 
position state consists of k zeroes and {N — k) ones, 
and the second term is obtained by replacing every zero 
and one of the first term with one and zero, respec- 
tively. The index k e {l,N/2} (iV=even) or fc e 
{1, {N- l)/2} (iV=odd). If N is even, each k group has 
2(^) members, except when k = N/2, for which there 

are (^2) members. If N is odd, there are always 2(^) 2Afe,jv-/c = | |*fc ) + | ) = "fcj + 



an 



(12) 



and 



a 



members in every fc-group. 

The final density matrix is diagonal in the cat ba- 
sis. This can be easily seen from the action of the noisy 
channel on the qubits. For example, a maximally entan- 
gled state of two qubits after the action of this channel 
becomes a mixed state which is diagonal in the Bell basis. 
The final density matrix can be written in the following 
form: 



= a% 



k,N-k 



k„N-k 



kj 

(13) 



where j^*^) are the members of the fc group. 

We can now state our basic result, which follows di- 
rectly from Eq. ((1(11) . 

Proposition 1 The state p^ is entangled if there is at 
least one k such that 



kj 



kj 



A > 2Afe 



k,N-k- 



(14) 



Suppose that there is a fc = fco such that the above 
inequality is satisfied. This means that the state p^y is 

(9) 

NPPT across the bipartite partition where fco parties are 
on one side and [N — fco) on the other [which we denote 



4 



by ko : (N — ko)]. Let us also note that choosing a specific 
j (i.e., a given permutation) implies choosing a specific 
set of k parties. As A is independent of j, the entan- 
glement condition does not depend on which parties we 
have chosen for a given number of parties k on one side 
of the bipartite cut. 

We now prove a useful property of the final state: 

Lemma 1 // the state pj^ is entangled for some bipartite 
partition k : {N — k), then it is entangled for all bipartite 
partitions m : (N — m) when m > k and 1 < fc, m < 
N/2{Neven),{N ~ l)/2{Nodd) . 

Proof. Let us first note that we are considering only 
bipartite partitions where k{m) parties are on one side 
and N — k{N — m) are on the other. This automatically 
puts the second constraint on the allowed values of k, m. 
We need to prove that \k,N-k > ^m,N-m, i.e., a''b^~'' + 

Ijk^N-k ^ (fnijN^m _^ ^m^W-m^ ^^^^ m > fc and 1 < 

fc, m < N/2{Neven), {N - l)/2, N = odd. Without any 
loss of generality we can assume that a > b because both 
sides are equal when a=b. Suppose, the inequality is 
not valid, i.e., a'^b^-'' + b^a'^''' < a™6^-" + 6"a^-™. 
Then factoring a common a'^b'^ from both sides we have 

k + l. Let f =a;. Then x^-^fc ^ ;l ^^i^^N-2k-i ^ ^^ich 
implies x^^"^^ _^N-2k-i < 2;' — 1. This can be rewritten 
as x^~'^^~^{x^ ~ 1) ^ a;' — 1 which is a contradiction 
because x > \ and N — 2k — l = N— {k + m) > (since, 
k + m<N). ■ 

It follows immediately from this lemma that: 

Corollary ^ (i) If the final state is entangled across the 
\ : N —1 cut then it is entangled across all other bipartite 
cuts, (ii) The most robust bipartite cut is N/2 : N/2 (N 
even) or {N±l)/2 : (TV^ l)/2 (N odd). 

Let us now consider an M-partition with the parti- 
tions labeled as d, G2, .., Gm, 2 < M < TV. We call 
such a choice an M-partition configuration. Let \Gk\ 
be the number of particles in the group Gk. Then an 
Af-qubit cat state can be distilled iff the state is NPPT 
for all possible bipartite cuts in the A/-partition config- 
uration. That includes the bipartite partitions \Gk\ ■ 
{N — \Gk\) VGfc, as well as the bipartite partitions ob- 
tained by combining a subset of the partitions Gm- Let 
|Gi| = miufe {|Gfe|}. Then we have the following lemma, 
that states under which condition one can distill an Af- 
qubit cat state. 

Lemma 2 An M-qubit cat state can be distilled if and 
only if A> 2A|G^|^Ar„|G,| . 

Proof. Using the same calculation as in Lemma ^ note 

that if A > 2A|G.|_Ar-|Gih ^ > '^^Gkl^N-iOkl ^^k- 
This guarantees that the final state is NPPT across all 
other possible bipartite cuts in the Af-party configura- 
tion. Now note that our state can be brought to a depo- 
larized form by local operations while satisfying the con- 
dition A > 2A|G.|_jv_|Gi| as well as preserving the NPPT 



property. For such a depolarized state distillability is 
guaranteed since the satisfaction of the above condition 
is both necessary and sufficient for Af-qubit distillability 
Hence our state is also distillable. ■ 
Note that the symmetry properties of the state and 
the hierarchical property of the A's guarantee that it is 
sufficient to compute the partial transposition of only 
one bipartite cut. The next lemma shows the connection 
between a bipartite cut and the largest M one can have 
in Af-qubit distillability. Let Q{N, fc) =quotient of (TV — 
fc)/fc. 

Lemma 3 For any given k in a bipartite cut fc : (A^ — fc), 
if the distillability condition A > 2Xk^N~k is satisfied, 
then the density matrix p'j^ is M-qubit cat state distil- 
lable, where maxAf = N/k if N is divisible by fc or 
maxAf = 1 + Q(-/V, fc) if N is not divisible by fc. 

Proof. Let fc be the smallest integer such that A > 
2Afcjv-fc. Clearly if N is divisible by fc, the maximum 
number of groups one can have is N/k. However if 
N is not divisible by fc, then the remainder of N/k is 
less than fc and hence 1 -I- N/k cannot be the number 
of optimal groups, as the partition {remainder of N/k) 
■.{N — remainder of N/k) is not NPPT. Therefore the op- 
timal number of partition must be 1 -I- Q{N, fc). ■ 

It is clear that every bipartite cut contains distillability 
information about some M-partition. In what follows we 
will analyze the entanglement threshold conditions for 
bipartite cuts. The previous lemmas guarantee that the 
bipartite threshold conditions are sufficient to describe 
distillability of the final state p^. 

VI. ENTANGLEMENT PROPERTIES OF 

In this section we study in detail entanglement proper- 
ties of the state /djy- noted before we will concentrate 
on the conditions for the state to remain entangled across 
the bipartite cuts. Let us first point out the effect parity 
(in the sense of odd/even number of qubits) can have on 
the inequality. As mentioned before, the parameters c, d 
can be either positive or negative. Without loss of gen- 
erality suppose |c| > \d\. When c, d have the same sign 
we have [recall Eq. (|12|l ] A = |c|^ -I- jdj^ irrespective of 
N being even or odd. However if c, d have opposite signs 
then A — \c\^ -f \d\^ when A^ is even or A = \c\^ — \d\^ 
when N is odd. This suggests that parity can have a dra- 
matic effect on the sufficient condition for entanglement. 

A. Sufficient condition for entanglement when c, d 
liave tlie same sign 

1. Exact condition for finite N 

Let us rewrite our basic NPPT condition H14|l explicitly 
as (|c|^ + \df) > a%^-^ + b'^a^-^. Let us note first 
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that, since from Eq. (O it follows that a > c and b > d, 
the NPPT condition is violated when a = b. 

From hereon we assume a ^ b. Note that if |c| = |d| 
the entanglement condition can still be satisfied if a 7^ 
b. Letting k — Na, where i < a < i, one obtains 



icra 



') > 5"^a(i-")^(l + (i)^'-'"^"^), where 
one can assume without loss of generality that a > b and 
|c| > Taking the logarithm of both sides and dividing 
by N, one obtains 



log|c| 



TV 



log(l + 



AT 

) > log (b^a^^- 



N 



log(l 



(l-2Q)Ar 



)• 

(15) 



This inequality is an exact sufficient condition for entan- 
glement corresponding to a bipartite cut k : {N — k) . 



2. Asymptotic condition when the subsystems and system 
sizes become macroscopic 

Consider first the case when the subsystem size k be- 
comes macroscopic in the limit N ^ 00. This implies 
a remaining a constant as iV — > 00. The asymptotic 
condition can easily be obtained from Eq. H15|) : 



> 6", 



(16) 



Note that d dropped out in the asymptotic limit. We 
can rewrite this as a lower bound on a in the asymptotic 
limit, using the normalization condition a + b = 1: 



a > 



log iik) 



Remarkably, the asymptotic condition depends only on 
two parameters, a and c. Noting that a = 1/M, Eq. H17|) 
puts an upper bound on M. 

Since < a < ^ we find from Eq. 11()|) a particularly 
simple form of the entanglement condition in the case of 
the most robust partition, i.e. 



the case a — ^: 



> ab 



(18) 



We note here that an identical inequality can also be 
obtained by putting the constraint /(a, |c|) < 1/2. 



3. Robustness and size of the subsystems 

An interesting feature of Eq. I|15(l is that we can sepa- 
rate out the asymptotic part from the A'^-dependent term. 



Let us rewrite Eq. (|15|l as 

log|c| > log (6"a(i-") 

log(l - 



+ 



1 

TV 



{l-2a)N 



log(l 



d 


N ' 

) 


c 





(19) 



In this form it is easy to see that the asymptotic condi- 
tion is achieved even for finite N if and only if log(l + 

(i)^^-^"^")=iog(i+ifr 

is satisfied: 



), i.e., the following condition 




(20) 



For instance, when N is even, for a choice of |c| — \d\, one 
obtains a — 1/2. This shows, rather surprisingly, that 
the asymptotic condition is the same as that of any finite 
N if certain conditions are met. Consistency demands 
that, in such a case, Eqs. H17|l . (|2()|l should be satisfied 
simultaneously. This leads to the condition |c||d| > a6, 
which cannot be satisfied in general, except for the case 
when |c| = \d\. We therefore conclude that latter surpris- 
ing property holds only for equi-grouped partition, and 
only when N is even, because a ^ 1/2 if is odd. This 
therefore shows another effect of parity of the number of 
qubits. 

We now examine Eq. (|19|l more closely. Let us denote 



= log(l ^ 
Eq. lO as 



J. (1-2q)7V 



(I) 



) — log(l + f ) and rewrite 



log |c| > log (6"a(i- 



N' 



(21) 



(17) Let us first observe that ^ > if a > ^1 



log(f)^ 

Therefore if Eq. H21|l is satisfied for some finite A^ when 
fi is positive, then the asymptotic condition of entangle- 
ment (|17ll is automatically satisfied as well. As before, 
demanding the consistency requirement one can show 

' /(a, |c|) if \c\\d\ < ab, which, in 



that 



1 



< 



i°g(f). 

general, is always satisfied. We can therefore conclude 
that existence of M-group entanglement for some finite 
N automatically implies entanglement also in the asymp- 
totic limit when the size of the partitions become macro- 
scopic, keeping the number of partitions M a constant. 
On the other hand if the inequality H21|l is not satisfied 
for some A^, it might still be satisfied for some large A^, 
as the condition itself gets relaxed as A^ — > 00. We sum- 
marize these considerations as follows: 

Proposition 2 Multi-partite entangled states with fixed 
number of partitions M are more robust the larger is the 
dimension of the constituent subsystems. 
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4- Fixed number of members in a group while number of 
groups and the size of the system become large 

We now discuss the second scenario where we ahow the 
number of particles in a group to remain constant, i.e., 
k is fixed while both M, N are allowed to become arbi- 
trarily large. For simplicity we assume that each group 
contains the same number of qubits. We first rewrite 
Eq. (O as 



log|c| > — log 6+ (1 



f 

TV 



log(l 



) log a 

n' ^ 



-a 



) - log(l + 



N 



(22) 



Suppose the above inequality is satisfied for some choice 
of fc, N. However, in the limit iV — > oo the inequality 
reduces to the condition log|c| > log a, which is false. 
Thus: 

Proposition 3 If we allow both the number of partitions 
and the number of qubits to become large, inequality (22) 
ceases to be satisfied. 



B. Condition for entanglement wlien c, d have 
opposite signs 

Let us now consider the case of odd number of qubits 
and c,d having opposite signs, whence A — \c\^ — ■ 
First, note that the condition (|14l) is considerably tighter 
than before. Second, there is no distillable entanglement 
in this case when c = d as opposed to the same-sign case, 
irrespective of whether a, b are equal or not. 

Let us write the condition (|14|l in this case as: 



log|c| > log (6"a(i-" 
log(l - 



1 

iV 



(l-2Q)Af 



) - log(l 



d 


N ' 

) 


c 





(23) 



In the asymptotic limit, we obtain the same condition as 
before. However, notice that the term in the parenthe- 
ses on the right side of the inequality is always positive. 
This implies that as we increase the number of particles 



robustness always increases for any choice of a. This 
means that if the state is entangled for some choice of 
M,N, it always remains entangled in the large N limit, 
as robustness increases with N as long as the number of 
partitions M remain unchanged. 



VII. CONCLUSIONS AND OPEN PROBLEMS 

To summarize, in this work we have studied the ro- 
bustness of A^-qubit cat states under a rather general 
decoherence model. In this model every qubit is indepen- 
dently coupled to the environment. The noisy channel is 
described by a completely positive map with arbitrary 
probabilities assigned to the various errors described by 
the Pauli operators. Our findings show that macroscopic 
entanglement is more robust in higher dimensional sys- 
tems while keeping the number of spatially separated par- 
ties constant. We have also shown that states with even 
or odd numbers of qubits can have qualitatively different 
properties that are not observed in simpler noisy chan- 
nels such as depolarizing channel. Furthermore, we have 
shown that in the asymptotic limit the entanglement con- 
dition depends only on two noise parameters even though 
the noise model itself is described by three independent 
parameters. 

The present work focuses on the GHZ class of states. It 
would be interesting to see if the above observations hold 
true for other classes of multi-qubit pure states. Kempe 
and Simon have studied the robustness of W-states 
for three qubits, and other inequivalent classes of four- 
qubit entanglement, using a depolarizing channel. Their 
work has shown that GHZ states are more robust than 
the other classes of states. It would be interesting to 
test whether this holds true under the more general noise 
model considered here. The same comment applies to the 
class of graph states studied by Diir and Briegel 

Note added: Upon completion of this work we came to 
know about related work by Hein et al. |20| |. which uses 
a Markovian master equation approach to model deco- 
herence, whereas in this work we have used the (formally 
exact) Kraus operator sum representation. 
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